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We propose helium-4 spallation processes induced by long-lived stau in supersymmetric standard 
models, and investigate an impact of the processes on light elements abundances. We show that, 
as long as the phase space of helium-4 spallation processes is open, they are more important than 
stau-catalyzed fusion and hence constrain the stau property. 



I. INTRODUCTION 

Quests for the physics beyond the Standard Model 
(SM) will reach a new stage at the TeV scale. Among 
the expected interesting signals of the new physics are 
those provided by exotic charged particles (charged mas- 
sive particles; CHAMPs) with a long lifetime. The pres- 
ence of such particles is predicted in many notable mod- 
els beyond the SM, although its identity depends on the 
models one assumes. CHAMP hunting is indeed one of 
the major issues of the high energy experiments, and its 
collider phenomenology is enthusiastically studied [i"14l4j] ; 
it also motivates other researches including neutrino tele- 
scope observations and cosmology fl8i - l2lj |. 

Long-lived CHAMPs will play interesting roles in the 
Big-Bang Nucleosynthesis (BBN) as well. The light nu- 
clei will interact not only with the CHAMPs during the 
BBN processes [221 - I39I ] , but also with the decay products 
of the CHAMPs in the post-BBN era jijS3. The stan- 
dard scenario of the BBN will thus be altered, and so is 
the abundance of the light elements at the present time. 
One can thus constrain the models beyond the Standard 
Model by evaluating their prediction on the light ele- 
ments abundance and comparing it with the current ob- 
servations. We can then give stringent predictions for the 
forthcoming experiments and observations according to 
these constraints. 

The Standard Model extended with supersymmetry 
(SUSY) is one of the models that can accommodate such 
long-lived CHAMPs. With the i?-parity conservation, 
the lightest SUSY particle (LSP) is stable and become 
a cold dark matter. Interestingly, it can offer a long- 
lived CHAMPs if the LSP is the bino-like neutralino Xi- 
Coannihilation mechanism is required to account for the 
dark matter abundance in this case [48J, where the LSP 
and the next-lightest SUSY particle (NLSP) are almost 
degenerate in mass. Staus, denoted by f and a possible 



candidate of the NLSP, can acquire a long lifetime when 
the mass difference with the LSP is less than the mass of 
tau leptons. This is due to the phase space suppression of 
the final state that necessarily consists of three particles 
or more. Noting that such long-lived staus will be copious 
during the BBN [H, [Hoj], we have shown in [H, H2, S3] 
that their presence indeed alters the prediction of the 
standard BBN and possibly solve the discrepancy of the 
lithium abundance in the Universe through the internal 
conversion reactions. 

In this article, we improve our analyses by including 
new reactions of 
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in which (f 4 He) represents a bound state of a stau and 
4 He nucleus. Reaction ((T|) is essentially a spallation of 
the He nucleus, producing a triton t, a deuteron d, and 
neutrons n. Presence of such spallation processes has 
been ignored so far due to the naive expectation that the 
rate of the stau-catalyzed fusion [22j 



(f 4 He) + d^ f - 



'Li 



(2) 



is larger than the reaction (JXJ) . Indeed, the cross section 
of Eq. (JU) is much larger than that of 4 He + d -> 6 Li + 
7 by (6 - 7) orders of magnitude (5lj . 

We point out that this expectation is indeed naive; the 
reaction Eq. (fTJ) is more effective than Eq. @ as long as 
the spallation processes are kinematically allowed. The 
former reaction rapidly occurs due to the large overlap 
of their wave functions in a bound state. On the other 
hand, the latter proceeds slowly since it requires an ex- 
ternal deuteron which is sparse at the BBN era. The 
overproduction of t and d is more problematic than that 
of Li. This puts new constraints on the parameters of 
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overproduction of the various light elements due to these 
processes. 

In this section, we calculate the rate of the spallation of 
the 4 He nucleus. We compare the result with the rate of 
the stau-catalyzed fusion to show that the former is larger 
than the latter in a large part of the parameter space, and 
thereby show that the spallation plays a significant role 
in the BBN. 

The 4 He spallation processes of Eq. ([1]) is described by 
the Lagrangian 





C = rx° 1 (g L P L +g R P R )T 
+ V2G f v t ^P l tJ^ + h. 



(3) 



where G-p = 1.166 x 10~ 5 GeV~ 2 is the Fermi coupling 
constant, ^l(R) represents the chiral projection opera- 
tor, and is the weak current. The effective coupling 
constants and <?r are given by 
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where g is the SU{2)\ J gauge coupling constant and 6-w is 
the Weinberg angle. The mass eigenstate of staus is given 
by the linear combination of tl and tr, the superpartners 
of left-handed and right-handed tau leptons, as 



f = cos# t tl + sin# T e 17t tr. 



(5) 



FIG. 1: He spallation processes. 



Here 6 T is the left-right mixing angle of staus and "f r is 
the CP violating phase. 



the minimal supersymmetric standard model (MSSM). 
Note that there is no reaction corresponding to Eq. (fT]) 
in the gravitino LSP scenario [52j . 

The purpose of this work is to understand the im- 
pact of 4 He spallation processes on light element abun- 
dances. In Section HH we analytically calculate its re- 
action rates, and compare its timescale with that of the 
reaction Eq. @ . In Section IIII1 we calculate all of light 
element abundances including exotic reactions, i.e., the 
4 He spallation processes, the stau-catalyzed fusion, and 
the internal conversion processes. We show the MSSM 
parameter space in which we can reproduce the observed 
abundances of both dark matter and light elements in- 
cluding 7 Li and 6 Li. Section ITVl is devoted to a summary. 



II. SPALLATION OF HELIUM 4 

Two types of reactions are possible for the bound state 
of a stau and a 4 He nucleus: (1) the stau-catalyzed fusion 
and (2) the spallation of the 4 He nucleus. The property 
of stau is stringently constrained in order to evade the 



A. (f 4 He) ^xi' + ^r+t + n 

First we consider the process of Eq. (flaj) . The rate of 
this process is expressed as 



T tn |l/>| 2 • CTUtn' 

where \ip\ 2 stands for the overlap of the wave functions of 
the stau and the 4 He nucleus. We estimate the overlap 
by 

2= (Zam H e) 3 ; (?) 

where Z and toho represent the atomic number and the 
mass of He, respectively, and a is the fine structure con- 
stant. We assumed that the stau is pointlike particle and 
is much heavier than 4 He nucleus so that the reduced 
mass of the bound state is equal to the mass of 4 He nu- 
cleus itself. The cross section of the elementary process 
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for this reaction is denoted by av tn and calculated as 
av tn = crw((f 4 He) 



Xi^rtn) 
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Here Pi and Ei are the momentum and the energy of the 
particle species i, respectively. 

We briefly show the calculation of the amplitude of this 
process, leaving the full calculation in Appendix. The 
amplitude is deconstructed as 



M((f 4 Re 

= (tnx?i/ T |£ 
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(tn|JTHe) (X>r\j,\r 



Here we omitted the delta function for the momentum 
conservation and the spatial integral. The weak current 
consists of a vector current and an axial vector 
current as = + gxA^, where g A is the axial 
coupling constant. The relevant components of the cur- 
rents in this reaction are V° and A 1 (i — 1,2,3). We take 
these operators as a sum of a single-nucleon operators as 



V 



A^Y^r^W' 



qr a 



(10) 



where q is the momentum carried by the current, r a is the 
spatial coordinate of the a-th nucleon (a G {1,2,3,4}), 
and t~ and a % a denote the isospin ladder operator and 
the spin operator of the a-th nucleon, respectively. Each 
component leads to a part of hadronic matrix element: 



(tn|^°| 4 He) = \/2X tn , 
(tn| 5A A+| 4 He) = V2g A Mtn, 
(tn| 3A A-| 4 He) = -y/2g A M t 
(tn| 5A A 3 | 4 He) = -V2g A M t 



(11) 



where A^ = (A 1 ± iA 2 )/y/2. Given the relevant wave 
functions of a He nucleus, a triton, and a neutron in 
Appendix, we obtain the hadronic matrix element as 



M t 
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3 (dHc + a t ) 4 
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(<?t + q n ) 



3aHc 6(anc 
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(12) 



Here qt and q n are three-momenta of the triton and the 
neutron, respectively, and an e and at are related to the 
mean square matter radius i? ma t by 
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(13) 



TABLE I: Input values of the matter radius 7? ma t for d, t, 
and 4 He, the magnetic radius -R mag for p and n, nucleus mass 
mx, excess energy Ax for the nucleus X, and each reference. 



nucleus i? mat ( mag) [fm]/[GeV _1 ] 


m x [GeV] 


Ax [GeV] 






p 0.876 / 4.439 [53] 


0.9383 [58] 


6.778 x 10" 


-3 


[59] 


n 0.873 / 4.424 [54] 


0.9396 [58] 


8.071 x 10" 


-3 


[59] 


d 1.966 / 9.962 [55] 


1.876 [59J 


1.314 x 10" 


-2 


[59] 


t 1.928 / 9.770 [56] 


2.809 [59J 


1.495 x 10" 


-2 


[59] 


4 He 1.49 / 7.55 [5JJ 


3.728 [59] 


2.425 x 10" 


-3 


[59] 



We list in Table fl] input values of the matter radius for 
the numerical calculation in this article. 



The remaining part is straightforwardly calculated to 



be 



Kx? "r\m\ 2 = \{x>r\m? = 4G|| 5R | 2 ^V^. 



n n o n JTl^aEy ( -rf 

|( X V r | J± |f)| 2 =4G 2 M 2 ^ r h T ^ 



(14) 



where E v and pi, are the energy and the z-componcnt 
of the momentum of the tau neutrino, respectively. We 
assumed that the stau and the neutralino are non- 
relativistic. This equation includes not only all the cou- 
plings such as Gf, <?l, and (?r, but also the effect of the 
virtual tau propagation in the Fig. Q] Note here that <?l 
coupling does not contribute. This is because the virtual 
tau ought to be left-handed at the weak current, and it 
flips its chirality during the propagation since the trans- 
ferred momentum is much less than its mass. 

Combining hadronic part with the other part, we ob- 
tain the squared amplitude as 

|X((f 4 He)^x?Mn)| 2 

<15) 



Integrating on the phase space of the final states, we 
obtain the cross section as 



8 / 32 \ 3 ^ 2 
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Here A tn , fet, and fc n are denned as 

A tn = Sm + A Hc - A t - A n - E h , 



k t = y/2m t A t 



(18) 



fc n = v /2m n A tn , 
nerj 

is the binding energy of (f 4 He) system. 



B. (f 4 He)^x? + "r + d + n+n 



The rate of another spallation process of Eq. (jlbp is 
where Ax is the excess energy of the nucleus X, and E\> similarly calculated. The cross section is calculated to 

be 



192 



0"Udnn = — 7-9 tan #wsin e T G F A dnn 



2 4 m n m d 



2a d 



m f m 2 \a He (ad + flHe)' 



3/2 



*dnn ■ 



(19) 



where 
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Here Adnn, k n and Ai(i = 1 - 5) are defined as follows: 
Adnn = <5rn + A Hc - A d - 2 A n - Eb, 

k n = \/2TO„Adnn , 



The rate is then obtained in the same manner as Eq. ^ . 
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C. (f 4 He)^x? + ^+P + n + n + n 

The cross section of spallation process of Eq. (fTc|) is 
calculated to be 

8 /32tt 3 ^ 3/2 



o-fpnnn = ( ) g 2 tan 2 W sin 2 T (1 + 3^)G| 

7r J \ a 

x a Ho A^ r 
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where 
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where A pnnn and &n are defined as follows: 

A pnmi = 5m + A Hc - A p - 3A n - 

, (24) 

fc N = V 2m N A pIlnll . 

In this calculation, we assumed proton and neutron have 
an identical kinetic energy, and then the factor k p and fc n , 
which are introduced to factorize their kinetic energies, 
are also identical. &n is the identical factor, and here we 
took mjM = m n . 

The reaction rate is obtained in the same manner as 
Eq. @. 

D. Comparing the rate of spallation reaction with 
that of stau-catalyzed fusion 

We compare the rate of the spallation and that of the 
stau-catalyzed fusion. We first note that the rate of 
stau-catalyzed fusion strongly depends on the temper- 
ature [51], and we fix the reference temperature to be 
30keV. Staus begin to form a bound state with 4 He at 
this temperature, which corresponds to cosmic time of 
10 3 s. Thus the bound state is formed when the lifetime 
of staus is longer than 10 3 s. 

Figure [5] shows the timescale of the spallation processes 
as a function of 5m. The lifetime of free stau is plotted 
by a solid line. We took the reference values of nif — 
350GeV, sin0 r = 0.8, and 7 r = 0. The inverted rate of 
the stau-catalyzed fusion at the temperature of 30keV is 
also shown by the horizontal dashed line. Once a bound 
state is formed, as long as the phase space of spallation 
processes are open sufficiently that is 5m > 0.026GeV, 
those processes dominate over other processes. There f 
property is constrained to evade the over-production of d 




0.05 0.1 0.15 0.2 

Mass difference (GeV) 

FIG. 2: Timescale of spallation processes as a function of 
5m and the stau-catalyzed fusion at the universe temperature 
T = 30keV [55|- The lifetime of free f (solid line) is also 
depicted. Here we took mf = 350GeV, sin# T = 0.8, and 

7r =0. 



and/or t. For 5m < 0.026GeV, the dominant process of 
(f 4 He) is stau-catalyzed fusion, since the free f lifetime is 
longer than the timescale of stau-catalyzed fusion. Thus 
light gray region is forbidden due to the over-production 
of 6 Li. 

This interpretation of Fig. 2 is not much altered by 
varying the parameters relevant with f. First cross sec- 
tions of spallation processes are inversely proportional to 
TOf , and then the timescale of each process linearly in- 
creases as Wf increases. Thus, even when mf is larger 
than mf = 350GeV by up to a factor of ten, the region 
of 6 Li over-production scarcely changes. Next we point 
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out that our result depend only mildly on the left-right 
mixing of the stau. Indeed, cross section of the 4 He spal- 
lation is proportional to sin 2 6 T . Its order of magnitude 
will not change as long as the right-handed component is 
significant. 



III. LIGHT ELEMENTS ABUNDANCES AND 
ALLOWED PARAMETER SPACE 

We numerically calculate the primordial abundances of 
light elements including 4 He spallation processes and f 
catalyzed nuclear fusion. Then we can search for allowed 
regions of the parameter space to fit observational light 
element abundances. 

So far it has been reported that there is a discrepancy 
between the theoretical value of 7 Li abundance predicted 
in the standard BBN (SBBN) and the observational one. 
This is called 7 Li problem. SBBN predicts the 7 Li to H 
ratio to be Log 10 ( 7 Li/H) = —9.35 ± 0.06 when we adopt 
a recent value of baryon to photon ratio t) = (6.225 ± 
0.170]x 10" 10 (68% C.L.) reported by the WMAP satel- 
lite [60| , and experimental data of the rate for the 7 Li or 
7 Be production through 3 He + 4 He -5> 7 Be + 7 61] ( 7 Li 
is produced from 7 Be by its electron capture, 7 Be + e~ 
— > 7 Li + v e at a later epoch). On the other hand, the 
primordial 7 Li abundance is observed in metal-poor halo 
stars as absorption lines (62j. Recent observationally- 
inferred value of the primordial 7 Li to hydrogen ratio is 
Log 10 ( 7 Li/H) = — 9.63 ± 0.06 (63l for a high value, and 
Log 10 ( 7 Li/H) = -9.90 ± 0.09 [6| for a low value. (See 
also Refs. [65l - r67| for another values.) Therefore there 
is a discrepancy at more than three sigma between the- 
oretical and observational values even when we adopt 
the high value of [63j . This discrepancy can be hardly 
attributed to the correction of the cross section of nu- 
clear reaction [68, 69]. Even if we consider nonstandard 
astrophysical models such as those including diffusion ef- 
fects [TOj, l71|j_h might be difficult to fit all of the data 
consistently [72| ]. 

In Figs. [3] and [4] we plot the allowed parameter re- 
gions which are obtained by comparing the theoretical 
values to observational ones for the high and low 7 Li/H, 
respectively. Vertical axis is the yield value of f at the 
time of the formation of the bound states with nuclei, 
Yf = riffs (s is the entropy density), and horizontal axis 
is the mass difference of f and Xi- We have adopted fol- 
lowing another observational constraints on the light ele- 
ment abundances: an upper bound on the 6 Li to 7 Li ra- 
tio, 6 Li/ 7 Li < 0.046+0.022 [HI, the deuteron to hydrogen 
ratio, D/H=(2.80±0.20) x 10" 5 [z3, and an upper bound 
on the 3 He to deuteron ratio, 3 He/D < 0.87 + 0.27 [74|]. 

The solid line (orange line) denotes a theoretical 
value of the thermal relic abundance for staus [37j 
while keeping observationally-allowed dark matter den- 
sity n DM h 2 = 0.11 ± 0.01 (2 a) [13 as total abun- 
dance. For reference, we also plot the observationally- 
allowed dark matter density in the figures by a horizontal 
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FIG. 3: Allowed regions from observational light element 
abundances at 2 a. Here we have adopted the higher value 
of the observational 7 Li/H in [63|] denoted by ( 7 L1/H)h, and 
have plotted both the 2<r (thin line) and 3<r (thick line ) only 
for 7 Li/H. The horizontal band means the observationally- 
allowed dark matter density. We have adopted mf = 
350 GeV, sin# T = 0.8, and j T = 0, respectively. 



band. 

At around 8m ^0.1 GeV, we find that 7 Li/H can be 
fitted to the observational value without conflicting with 
the other light element abundances. 1 As shown in Fig.[3l 
it should be impressive that the relic density is consistent 
with the allowed region at Yf = 2 x 10~ 13 at 3 a in case 
of the high value of 7 Li/H in (63|. 



IV. SUMMARY 

We calculated primordial abundances of all of light el- 
ements involving the helium-4 spallation processes, the 
catalyzed fusion, and the internal conversion processes. 
Newly included in the present work is the spallation of 
the 4 He in the stau- 4 He bound state given in Eq.(l). This 
process is only present in the model which predicts the 
long-lived charged particles due to the phase space sup- 
pression with the weakly interacting daughter particle. 

We calculated the rate of the helium-4 spallation pro- 
cesses analytically, and compared it with that of cat- 
alyzed fusion. We found that the spallation of 4 He nuclei 
dominate over the catalyzed fusion as long as the phase 



1 See also [26 
7 Li/H. 
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where J M is the hadronic current. This interaction allows 
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FIG. 4: Same as Fig. [3] but for the lower value of observa- 
tional 7 Li/H reported in ,'64], which is denoted by ( 7 L1/H)l. 



space of the spallation processes are open and hence the 
property of long lived stau is constrained from avoiding 
the overproduction of a deuteron and/or a triton. In 
spite of these new constraints, we found that the lithium 
discrepancy and the dark matter abundance can be si- 
multaneously solved in the parameter regions presented 
in Fig. 1 
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Appendix A: Amplitude of the Spallation Reaction 
of "He 



The interaction relevant to our scenario is 

Ant =T*^(. 9L F L + 9R P R )T 

+ V2G f (v t ^P l t)J^ + H.c. . 



(Al) 



the spallation of the nuclei such as 

(f 4 He) -> xS + v T + t + n, 
(f 4 He) -»• x? + v T + d + n + n, 
(f 4 He) ->xi + ^V+P + n + n- 



(A2a) 
(A2b) 
(A2c) 



to take place. In this section, we illustrate how we cal- 
culate the amplitude of the processes of Eqs. (|A2j) . 



'He ->■ tn 



(A3) 



The amplitude of the process (|A2aj) is given by 

X((f 4 He) -> xXtn) 

= <tn| jfHe) (x>t\jAt). 

in which we define the leptonic matrix element by 
(X>r\j,\r) 

= V2G F (x? v r \ [vt^Plt] [t( 9 lPr + <&a)x?r] \f > 

(A4) 

and the hadronic matrix element by (tn| J^| 4 He) . We 
separately calculate these two matrix elements. 



a. Leptonic matrix element 

The leptonic matrix element is directly calculated un- 
der the following simplifications: 

• Neutralino is treated as non-relativistic particle 
since its mass m^o is much larger than the mass 
of nuclei. 

• The momentum of a virtual tau is negligibly smaller 
than its mass due to the assumption that the stau 
and the neutralino are nearly degenerate with the 
mass difference of O(10 - 100)MeV. 

A straightforward calculation leads to 

KxVr|jo|f)| 2 = Kx>r|j 3 |r>| 2 = 4GiM 2 ^|^, 



\(x>r\j±\r)\ 2 = AG 2 F \gn\ 2 " 



mi 



(A5) 



where Ei,Pi, and mi individually stand for energy, four- 
momentum and mass of particles. 



b. Hadronic matrix element 

Calculation of the hadronic matrix element requires 
the explicit form of the hadronic current and the wave 
functions of the nuclei. 
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We need the wave functions of initial helium, final 
triton, and nucleon. Building up these wave functions 
requires special attention to the symmetry. The wave 
function consists of spatial, spin, and isospin parts, and 
should be antisymmetric under the exchange of the two 
nucleons. The spin and isospin of the nucleus dictates 
the spin and isospin part of the wave function. We then 
arrange the spatial part so that the total wave function 
be antisymmetric under the permutation of the nucle- 
ons. We model the spatial wave functions by Gaussian 
functions in terms of Jacobi coordinates. 

Let us make a wave function of 4 He by this prescrip- 
tion. The spin and isospin parts of the wave function is 
constructed according to S = and I = 0, and turns out 
to be 

| 4 Hc) = 



2y/Q 



i P npn)(imi> + imt) - imt) - itm)) 

+ ipnnpX-itm) - imt) + itm) + imt» 
+ inppn) (-itm) - imt) + itm) + imt)) ( A6 ) 
+ |npn P )(itm) + imt) - imt) - itm)) 
+ |ppnn)(-itm) + imt) + um) - imt)) 
+ |nnp P )(-itm) + imt) + um) - imt)) 

The above wave function is antisymmetric under the ex- 
change of two particles. Thus the spatial part ought to 
be symmetric, and is constructed as 

,3 



z 3 
7T 



f r 2 , 2 , 2 , 

}j-aHc r 1 +r 2 +r 3 + 
~(ri + r 2 + r 3 + r 4 ) 2 | 



x exp< 



(A7) 



where an c parameterizes the nuclear radius and related 
to mean square matter radius i? ma t by 
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1 



16 (flmat)L 



(A8) 



We applied the non-relativistic normalization for the nu- 
clear wave function, which is 

J d 3 rid 3 r 2 d 3 r 3 d 3 r 4 \ipa e (ri, r 2 , r 3 , r 4 )\ 2 

x S 3 (^(r 1 +r 2 + r 3 + r 4 )^J =1. (A9) 

The wave function (|A7[) is independent of the coordinate 
of the center of mass since the initial 4 He is taken to be 
stationary. The complete wave function of 4 He nucleus 
is a direct product of Eqs. (|A6[) and (|A7|) . 

The wave functions of other nuclei are similarly ob- 
tained. For the triton, the spin and isospin part of the 



wave function is formed according to S = 1/2,/ = 1/2, 
and I z = —1/2; it is given by 



|tt) = ^=[|pnn)(|m)-mt)) 

-|npn)(|m)-|ltt)) 
+ |rmp)(|m)-|m)) 

|t|) = i=[|pnn)(|m)-|m)) 

-|npn)(|m)-|tlt)) 
+ |rmp)(|m)-|m)) 

Its spatial wave function is 
A( r i,r 2 ,r 3 ) 

= k^J eXP V Qt ' 3 (ri +r2+r3 ' 



(AlOa) 



(AlOb) 



x exp< — a t 



1 



(r x + r 2 + r 3 ) 



(All) 



where q t is the center-of-mass momentum of the triton 
and a t is related to the mean square matter radius as 



1 



1 



2 {Rm&t)t 



(A12) 



The wave function of the final neutron is simply taken to 
be 



(A13) 



with trivial spin and isospin parts. 

Using these wave functions, we can calculate hadronic 
matrix elements of Eq. (| A3|) as 



(tn|J°| 4 He) - V2M tn , 
(tn|J+| 4 He) =V2g A M tD , 
(tn|J-| 4 He) =-V2g A M tn , 
(tn|J 3 | 4 Hc) =-V2g A M tn , 



(A14) 



where Mtn is defined as 



Mtn 



1287T anoflt 
3 (aHc + at) 



3/4 











jexp 


3an c _ 


— exp 


3an c 6(an c + at) 



(A15) 



c. Amplitude 



Combining the hadronic matrix elements with the lep- 
tonic matrix elements, we obtain the square amplitude 
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of Eq. 



|M((f 4 Hc)^x> T tn)| 2 

= |(tn|J°| 4 He)| 2 \(x>r\jo\f)\ 2 

+ |(tn|J-| 4 Hc)| 2 |(x> T |j+|f)| 2 

+ |(tn|J+| 4 He)| 2 |(x?^|i-|r)| 2 

+ |(tn|J 3 | 4 He)| 2 |<X?^|i3|f)| 2 

8m^oGp|gR| 2 
= Xl J (l + 3gi)Kn^ 



(A16) 



2. 4 He ->• dnn 

The calculation presented in the previous section is 
applicable to another spallation process of Eq. (|A2b[) . 
The amplitude we need is 



X( 4 Hc^dnn) 
= (x>r\j^f)(dnn\J^Re) 



(A17) 



The leptonic matrix element is same as in Eq. (|A3p and is 
already calculated in Eq. (|A5[) . while the hadronic matrix 
element requires a calculation anew. 



a. Hadromc matrix element 

We need the wave function of the final states, which in- 
clude a deuteron and two neutrons. The spin and isospin 
of a deuteron is S = 1 and / = 0, and the corresponding 
wave functions are 



| d) +l} = -^(|pt}|nt)-|nt)|pt» 
|d,0) = i(|pt)|n|)-|nt)| P |) 



|p|)|nt)-|n|)|pt)) 



f 



|d,-l) = - 7 =(|p|)|n|)-|n|)|p|)). 
The spatial part of the wave function is given by 
ipd{ri,r 2 ) = [-) 

-o d (ri - r 2 ) 



(A18) 

(A19) 
(A20) 



x exp I iqd • j exp 



(A21) 



where <7d is the center-of-mass momentum and ad is re- 
lated to the mean square matter radius as 



f 



«d 



8 (R 



mat ) d 



(A22) 



The spin of two neutrons can be S — and 5=1. For 
each case, spin and isospin parts of the wave functions 
are 



no 



(|nt>|n|) -|n|)|nt> 



V2 

|m,+l> = |nt)|nt), 

|m,0) = -L(|nt)|n|) + |n|)|nt>, 

|ni,-l> = |n|)|n|), 



(A23) 
(A24) 
(A25) 
(A26) 



where |nj) expresses spin and isospin part of wave func- 
tion of S = i. Since the spin and isospin part of S = 1 
are symmetric under the exchange of two particles, the 
spatial part of the wave function ought to be antisym- 
metric. The spin and isospin part of S — 0, on the other 
hand, is antisymmetric under the exchange of two parti- 
cles, then the spatial part of the wave function ought to 
be symmetric. Therefore the spatial parts of each wave 
function are given by 



V2 



{exp [i(q n i ■ n + q n2 ■ r 2 )] 
+ exp [i(q n 2 ■ n + q nl ■ r 2 )]} : 



(A27) 



^ni(ri,r 2 ) = ^={exp [i(q n i • n + q n2 ■ r 2 )] 



V2 



exp [i(q n2 ■ ri + q nl ■ r 2 )]} . 



(A28) 



Using these wave function, we can calculate Hadronic 
matrix elements of Eq. (|A17[) as 



(dm|J°| 4 He) 


= VSMdnn , 


(d ni |J+| 4 He) 


= -V2gAMdnn 


(dm|J-| 4 He) 


= \/2gxMdTm , 


(dm|J 3 | 4 Hc) 


= V2.9A.Mdnn , 


(dn |J+| 4 Hc) 


= V2 3A X d „n , 


(dn |J-| 4 He) 


= -V2.9A.Mdnn 


(dn |J 3 | 4 He) 





(A29) 
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where Mdnn and M A „„ are defined as 



M dl m 



327r 2 a d 



aHe(aHo + ad) 2 
1 / 327r 2 a d 



exp 



4gg 2 + 4g n2 ■ g d + 3g 2 

8aHo 



exp 



Mli + 4gni ■ g d + 3gg 

8anc 



\ a He(aHc + Od) 2 



exp 



4g 2 2 + 4g n2 ■ g d + 3g 2 

8ohc 



exp 



4ggi + 4gm ■ g d + 3gg 

8anc 



2 exp 



Sg 2 ! + 2g ni ■ g n2 + 3g 2 2 1 (g d + q n i + g n2 ) 



8aHc 



(A30) 



b. Amplitude 

Combining the hadronic matrix elements with the lep- 
tonic matrix elements, we obtain the square amplitude 
of Eq. (jXTTT) as 

|X((f 4 Hc) -> x?^rdnn)| 2 

= |(dnn|J°| 4 He)| 2 \{x>r\j \r)\ 2 

+ |(dnn|J-| 4 Hc)| 2 Kx?i/ r |i + |f)| a 

+ |(dnn|J+| 4 Hc)| 2 \(x>r\j-\r)? 

+ |(dnn|J 3 | 4 He)| 2 \{x>r\h\f)\ 2 

12myoG||g R | 2 on n 

= Xl 2 ((1 + 2gl)Ml nn + 2glM'ln)E v ■ 

(A31) 



where we set e^u = +1. Here indices i, j, k, I are corre- 
sponding to spacial coordinates, i.e., spacial part of the 
wave function for |n t)i| n t)il n -J-)fc|p i)l is 



exp 



«(g n i • n + g n2 • rj- + g n3 • r fc + g p • n) 



(A34) 



Here each spacial part of proton and neutron are simply 
taken to be plane wave form. Similarly, the spin and 
isospin parts of the system for S = 1 and S z = { — 1, 0, 
+1} are given by 



| pnnn (5 = 1, S z = +1)) 

= ^cy W |nt)<|nt)i|n4.)fc|pt>J 



(A35) 



3. 4 He — > pnnn 

The matrix element for the 4 He — > pnnn process is 

(A32) 



|pnnn(S f = 1,S Z = 0)) 



|n t)i|n t)i|n |) fe |p l)i - |n |>,|n |)> f) fc |p t)," 

(A36) 



M{ 4 Hc — > pnnn) 



= (XiVr|jMl^)(pnnn|^| 4 He) . 

The calculation is also performed with an identical step 
as that of other 4 He spallation processes. 



a. Hadronic matrix element 

The final state of the process is a system composed of a 
proton and three neutrons, and two types of the systems 
could be brought; (1) S = and S z — via vector current 
(2) S = 1 and S z = {—1, 0, +1} via axial vector current. 

The spin and isospin part of the system for S = and 
S z = is given by 

|pnnn(5 l = 0,S Z = 0)) = -^j=e ijk i 

|n t)i|n t» DfclP i)i + |n i) t \n |),|n f>fe|p t)i 



|pnnn(5 = 1,S Z = —1)) 

1 (A37) 
= ^y|^ife/l n l)>l)il n t)fc|p|); 

Spatial parts of them are same as the system for 5 = 
and S z — 0. 

The hadronic matrix element in (IA32[) is calculated 
with built wave functions and explicit form of each cur- 
rent as follows, 

(pnnn|J°| 4 He) = V2M pnnn , 
(pnnn|J + | 4 Hc> = -V2g A M P nnn, 
(pnnn|J+| 4 He) = V2g A A^ P nnn, 
(pnnn | J 3 | 4 He) = V2 ffA A^p„nn, 



(A33) 
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where VW pn nn is define as follow: 




(A39) 



r 



b. Amplitude 

Combining the hadronic matrix elements with the lep- 
tonic matrix elements, we obtain the square amplitude 
of Eq. (TM21) as 

|.M((f 4 He) -> x?^rpnnn)| 2 

-|(pnnn|J°| 4 Hc)| 2 |<x$"r|jo|f>| 2 

+ |(pnnn|J-| 4 He)| 2 |<xK|j+|f>| 2 

+ |<pnnn|J+| 4 He)| 2 \{x%\j-\f)\ 2 (A40) 

+ |(pnnn|J 3 | 4 He)| 2 \(x>r\h\f)\ 2 

8m^oGp|gR| 2 
= Xl J (l + 3g%)Ml nnn E u . 
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